First-principles study of the optical properties of Mga;Ti(i_a;)H2 



00 

o 
o 

(N 

u 

< 



C/3 



M.J. van Setten\ S. Er^, G. Brocks^, R.A. de Grooti-^ and G.A. de Wijs^ 

^Electronic Structure of Materials, Institute for Molecules and Materials, 

Faculty of Science, Radboud University Nijmegen, 

Toemooiveld 1, 6525 ED Nijmegen, The Netherlands 

^Computational Materials Science, Faculty of Science and Technology and MESA+ Institute for Nanotechnology, 

University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands and 

'^Laboratory of Solid State Chemistry, Zemike Institute for Advanced Materials, 

Rijksuniversiteit Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands 

(Dated: April 2, 2008) 

The optical and electronic properties of Mg-Ti hydrides are studied using first-principles density 
functional theory. Dielectric functions are calculated for Mga:Ti(i_^)H2 with compositions x = 0.5, 
0.75, and 0.875. The structure is that of fluorite TiH2 where both Mg and Ti atoms reside at the Ti 
positions of the lattice. In order to assess the effect of randomness in the Mg and Ti occupations we 
consider both highly ordered structures, modeled with simple unit cells of minimal size, and models 
of random alloys. These are simulated by super cells containing up to 64 formula units (Z = 64). 
All compositions and structural models turn out metallic, hence the dielectric functions contain 
interband and intraband free electron contributions. The former are calculated in the independent 
particle random phase approximation. The latter are modeled based upon the intraband plasma 
frequencies, which are also calculated from first-principles. Only for the models of the random alloys 
we obtain a black state, i.e. low reflection and transmission in the energy range from 1 to 6 eV. 



PACS numbers: 71.20.-b,71.15.Nc,61.72.Bb,74.62.Dh 
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I. INTRODUCTION 

Since the discovery of the switchable mirror YH^; by 
Huiberts et ali^ in 1996 several other metal hydride 
systems that behave as switchable mirrors have been 
discoveredjii^'^'^'^'^. The metals are reflective, but af- 
ter hydrogenation become semiconductors and hence in 
most cases become transparent. Especially when an alloy 
with high hydrogen mobility is used and when applied as 
thin films the optical switching can be fast, reversible and 
robust >^ 

Recently, meta-stable thin films composed of various 
ratios of magnesium and titanium have been shown to 
exhibit remarkable optical properties which could be es- 
pecially useful for smart solar cell coatings and hydrogen 
sensor applicationsi^i^ In the dehydrogenated state the 
films are highly reflective. Upon hydrogenation they be- 
come black, i.e. have a low reflection and high absorption 
in the energy range of the solar spectrum. The struc- 
tural and electronic characteristics of this black state are 
a topic of intensive research^^iiS. 

Obtaining structural data for these systems is difficult. 
Single crystals cannot be grown as the compounds are 
thermodynamically unstable. For a 7:1 ratio of magne- 
sium to titanium Kyoi et ai, using a sample prepared at 
high pressure, could determine the crystal structure.!^''' 
It is similar to the fluorite Till2 structure. Notten and 
co-workers (Refs. [lulls') and Borsa et al^^ have shown 
that for higher titanium content (x < 0.8) the structure 
is fluorite-like as well. At lower lower titanium concen- 
trations a rutile (a-MgH2 like) structure is found. Inter- 
estingly, the kinetics of the hydrogen ab/dcsorption re- 
actions are much faster in the fluorite structures than in 
the rutile structures'^ The equality of the molar volume 



of TiH2 and Mg has been used to explain the structural 
stability of the meta-stable "phases" M- Calculations using 
density functional theory (DPT) find the same composi- 
tion dependence of the relative stability of the rutile and 
fluorite structures f^^ 

The origin of the "black state" is not understood. 
Its explanation is, of course, intimately related to the 
electronic structure of the hydride. In experiments, 
usually less than two hydrogen atoms per metal atom 
can be reversibly absorbed and desorbedi^ Electrochem- 
ically 1.7 hydrogen atoms per metal atom can be stored 
reversiblyji^ Moreover, the crystal structure of Mg7TiHa; 



as determined by Kyoi et 



is estimated to contain 



1.6 hydrogen atoms per metal atom. All this suggests 
that maximally two metal electrons per metal atom can 
be transferred to hydrogeni^ Hence the Ti atoms re- 
main in an open-shell configuration with at least two 
d-electrons left. The above reasoning has been con- 
firmed by DFT calculations on crystalline Mga;Ti(i_2,)H2 
structures;^ The calculated densities of states show pre- 
dominant hydrogen-s character below the Fermi level, 
which is typical for metal hydrides, and titanium-d states 
at the Fermi level. These are likely to form metallic 
bands, so one expects Mg2:Tii_^H2 to be reflective in- 
stead of black. Moreover, experimental investigation of 
the electrical transport properties reveal high resistivity 
1.32-1.9 mricm, with a logarithmic, hence non-metallic, 
temperature dependence.^ In order to explain these re- 
sults the formation of a so-called coherent crystal struc- 
ture was proposed;^ wherein regions of insulating Mgll2 
and metallic TiH2 coexist. 

In an effort to advance the understanding of the black 
state, in this paper we report a computational study of 
the optical and electronic properties of Mg2;Ti(i„2,)Il2 for 



X = 0.5, 0.75 and 0.875. We employ both simple (i.e. min- 
imal) unit cells and large super cells of the same compo- 
sitions. The latter model random alloys wherein the Mg 
and Ti are distributed over the lattice sites of a TiH2- 
like structure. Details of these models are presented in 
Ref. M- 

The dielectric functions consist of intra and interband 
contributions: 



£{uj) = £intor(w) + £intra(w) 



(1) 



These are calculated separately. For the interband con- 
tribution Sinter wc use first-principlcs DFT in the inde- 
pendent particle random phase approximation. Since the 
materials in question are metals, the intraband dielectric 
part Eintra does not vanish. It is modeled from the plasma 
frequencies u)-p, which are calculated from first-principles 
as well. 

In Sec. HI] the technical details of the calculations are 
summarized. Sec. lIIII contains a concise description of the 
structural models used. Results on the dielectric func- 
tions are presented in Sec. IIVI Finally, in the discussion 
section (Sec. |V| results are compared with experiment.— 



II. COMPUTATIONAL METHODS 

First-principles DFT calculations were carried out with 
the Vienna Ah initio Simulation Package (VASP),— li^ii^ 
using the projector augmented wave (PAW) method. ^^J^ 
A plane wave basis set was used and periodic bound- 
ary conditions applied. The kinetic energy cutoff on 
the Kohn-Sham states was 312.5 eV. For the exchange- 
correlation functional we used the generalized gradi- 
ent approximation (GGA) ^ Non-linear core corrections 
were applied. ^-'^ 

The Brillouin zone integrations were performed using 
a Gaussian smearing method with a width of 0.1 eV.— . 
The k-point meshes were even grids containing F so that 
the band extrema are typically included in the calculation 
of the dielectric functions. The convergence of the dielec- 
tric functions and intraband plasma frequencies with re- 
spect to the k-point meshes was tested by increasing the 
number of k-points for each system separately. A typical 
mesh spacing of about 0.01 A~^ was needed to obtain 
converged results. 

The calculations of the complex interband dielectric 

functions, einter(w) = ei„ter('^) + ^Sinter(^)' "^'^^^ P'^'^" 

formed in the random phase independent particle approx- 
imation, i.e. taking into account only direct transitions 
from occupied to unoccupied Kohn-Sham orbitals. Local 
field effects were neglected. The imaginary part of the 
macroscopic dielectric function, £\J^^Juj), than has the 
form: 



e|^tcr(q.^) = 
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FIG. 1: (Color online) Dielectric functions of a free electron 
gas according to equation (3] 



where q denotes the direction of q and w, k -I- q and c, k 

label single particle states that are occupied and unoccu- 
pied in the ground state, respectively, e, u are the single 
particle energies and the translationally invariant parts 
of the wave functions. V is the volume of the unit cell. 
The real part, £i,itor('^)' ^^ obtained via a Kramers-Kronig 
transform. Further details can be found in Ref. 2^ 

Most optical data on hydrides are obtained from micro- 
and nano-crystalline samples whose crystallites have a 
significant spread in orientation. Therefore the most rele- 
vant quantity is the directionally averaged dielectric func- 



tion, i.e. e 



(2) 



(oj) averaged over q. 



The intraband dielectric function, £intra('^) — 
e[J^^^{uj) -f ieiJ^^^iijj), is calculated from the free electron 
plasma frequency LOp. 



e|,itra(^) 



= 1- 



^^ 



-(2) 



.(^) 






■ ix'7^ 



(3) 
(4) 



where an inverse lifetime, 7 — 0.01 eV, is used. For 
7 = the reflection would be perfect up to the plasma 
frequency and zero beyond. Finite values of 7 decrease 
the reflection below Wp and smoothen the reflection edge 
at Lo-p. For metals 7 values are in the order of 0.1 eV. By 
using a small value the influence of the interband part is 
emphasized. Calculating 7 from flrst principles obviously 
goes beyond DFT. For three values of 7 the free electron 
eintra is plotted in Fig. [Ij 

The plasma frequency uj-p is calculated as an integral 
over the Fermi surface according to: 



47re^ 
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(5) 
with 5k the weight factors, and /(cnk) the occupation 
function. Again we use directionally averaged values. 
Further details on the calculation of the plasma frequency 
can be found in Ref. 24. 

Finally the optical constants, the refraction index, n, 
and the extinction coefficient, k, and the absorption. A, 



reflection, R, and transmission, T, are calculated using 
the standard expressions: 



^\l ^''' + \fm" + i^^'^y 



K 

A 
R 
T 



1 — exp(—KUjd/c) 

(n-l)2+K2 

{n + 1)2 + k2 
(l-i?)(l-A), 



(e(2))^ 



(6) 

(7) 

(8) 

(9) 

(10) 



with e*-^^ and e^"^' the real and imaginary part of e, d the 
slab thickness and c the speed of light in vacuum. The 
reflection and transmission spectra are constructed to 
simulate the substrate/Mgj:Ti(i_3.)H2/palladium setup 
as was used in the experiments by Borsa et al^ All inter- 
nal reflections and absorptions in the three layer system 
are taken into account. 



III. STRUCTURES 

The calculation of the dielectric functions is performed 
using the crystal structures developed in Ref. J^. A brief 
summary of compositions and cell parameters is given in 
Table H] In short they were constructed in the following 
way. 

In the case of a; = 0.875 the simple cell is just the op- 
timized experimental cell with composition Mg28Ti4ll64. 
For X — 0.5 and 0.75, two and three atoms, respectively, 
out of the four titanium atoms in the conventional fee 
TiH2 cell were replaced by magnesium. Thus the unit 
cells have compositions Mg2Ti2H8 and MgaTiHg respec- 
tively. 

To simulate the random alloys super cells are used. 
These are also based on the fluorite Ti4ll8 [Z — \) build- 
ing block. For x = 0.5 and 0.75 2x2x2 super cells were 
constructed and for x = 0.875 a 2 x 1 x 1 super cell. Again 
Ti were replaced by Mg, but now such as to approximate 
random alloys most efficiently (see Ref. ilO for details) . 

For all the models constructed the positional and cell 
parameters were relaxed. The cells remain close to cubic, 
see Table m The angles between the crystal axes are close 
to 90°, except for the x = 0.75 and 0.875 super cells 
where there is a small deviation. 



IV. DIELECTRIC FUNCTIONS 



TABLE I: Number of formula units (Z), lattice parameters 
(A) and shortest Ti-Ti interatomic distances (A). The first 
row for each composition contains the data of the simple cell 
and the second of the super cell. 



a 



Ti-Ti 



0.5 

0.75 

0.875 



4 
32 
4 
32 
32 
64 



4.72 
9.08 
4.62 
9.29 
9.36 
18.75 



4.72 
9.03 
4.62 
9.30 
9.36 
9.42 



4.65 
9.09 
4.62 
9.27 
9.36 
9.73 



3.16 
3.15 
4.62 
3.13 
6.62 
3.08 
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FIG. 2: (Color online) Imaginary parts of the interband di- 
electric functions of Mga;Ti(i_^)H2 in the primitive cells. 
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FIG. 3: (Color online) Imaginary parts of the interband di- 
electric functions of 'iAgxT^{i-x)'^2 in the super cells. For 
X — 0.5 the dielectric function peaks to a value of 130. 



A. Interband dielectric functions 

Figures [2] and [3] show the calculated imaginary parts of 
the interband dielectric functions in the simple and su- 
per cells respectively. In general the dielectric functions 
exhibit a peak at low energy, below 2 eV, followed by a 



relatively flat tail. In the super cells the dielectric func- 
tions have higher peaks and flatter tails. An exception 
is the dielectric function of the x = 0.875 super cell. It 
does not have a peak at low energy. 

To understand the main features of the dielectric func- 
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FIG. 4: (Color online) Electronic densities of state of 
the simple and super cells for the three compositions of 
Mg3;Ti(i_^)H2. The zero of energy is taken at the Fermi level. 



tions we study the densities of states (DOS', Fig. \^. 
Below the Fermi level the DOS' have an approximately 
triangular shape with a predominant hydrogen-s and 
niagnesium-s and p character. At the Fermi level the 
DOS' are dominated by titanium-d states. Above the 
Fermi level the DOS' have a mixed character of hydrogen- 
s and p, magnesium-s and p and titanium-d. As an illus- 
tration the angular momentum projected partial DOS' of 
the simple cell with x = 0.75 is shown in Fig. [5] 

The electron density is localized at both the hydro- 
gen and titanium atoms. This is illustrated in Figure [SI 
where the electron localization function (Ref. ^5) is plot- 
ted, also for the simple cell Mgo.75Tio.25H2. At the tita- 
nium atoms the only available states are of d-character. 
Since d — d transitions are very weak the main contribu- 
tion to the dielectric function is from s—p transitions on 
the hydrogen atoms. 

In first order the imaginary part of the dielectric func- 
tion can be described as the joint density of states 
(JDOS) divided by w^. The hydrogen DOS has a dip near 
the Fermi energy. It increases when moving away from 
the Fermi energy, both to lower and higher energies. This 
causes the JDOS at the hydrogen atoms to increase more 
than linearly (in the region from to about 8 eV) . When 



FIG. 5: (Color online) Angular momentum projected 
partial densities of states (PDOS) [states/(atom eV)] of 
Mgo.75Tio.25ll2 in the simple unit cell. The zero of energy 
is at the Fermi level. The PDOS' are calculated in spheres 
with radii of 1.3, 1.3 and 0.8 A for Mg, Ti an H respectively. 




FIG. 6: (Color online) Electron localization function of 
Mgo.75Tio.25H2 in the simple unit cell. The Ti atoms are at 
the corners, the Mg atoms at the face centers and the H atom 
(small) are inside the cell. The electron localization function 
is plotted on a (110) cut-plane. Localization at the H and Ti 
atoms is evident. 



TABLE II; Squared plasma frequencies, 
band transitions. 



jp (eV^), from intra- 



simple cell 



super cell 



2; = 0.5 
X = 0.75 
X = 0.875 



10.7 
12.8 
8.6 



1.3 
0.4 
0.1 



divided by w^ this increase is rather effectively compen- 
sated. Hence the dielectric function does not vary much 
in the interval from 2 to 8 eV. This reasoning applies 
for all titanium concentrations that we considered. In 
the super cells the averaging over the various hydrogen 
DOS' gives rise to a further leveling out of the dielectric 
functions. Of course vibrational effects, that are lacking 
in our K calculation, will even further smoothen the 
dielectric functions. 

The peaks at the lower end of the energy range do 
arise from d-d transitions on the Ti atoms. Although 
their oscillator strength is small, division by w^ causes 
them to stand out nevertheless. For x = 0.875 the Ti 
DOS at i^Fcrmi is strougly suppressed in the super cell 
(see Fig. 2]). This correlates with the absence of the peak 

in ejntCT foi' the super cell at this composition (Figs. [H [3]). 
At the other two compositions, we see no clear correlation 
between difference in DOS and peak shape (comparing 
simple and super cells). The higher peaks in the DOS in 
the super cells can be understood when we realize that 
the effective "back folding" of the d bands and their mu- 
tual interactions (caused by the randomization) results 
in a flattening of the bands. Some of these bands will 
be very close to i^Pormi and their transitions will thus 
be "boosted" , both by the flatness of the bands and the 
small transition energy. This discussion anticipates the 
discussion in the next section. 



B. Intraband plasma frequencies 

The intraband plasma frequencies, which have been 
calculated according to Equation^ are listed in Table [TTl 
The squared plasma frequencies from the super cell cal- 
culations are between one and two orders of magnitude 
smaller than those of the simple cells. Therefore the 
edge on the free carrier reflectively occurs at consider- 
ably lower energies in the random alloys. In our models 
the maximum cUp is 1.1 eV. This is, however, for only 
one realization of a random model at x — 0.5 and it is 
well conceivable that calculations on large models could 
yield even lower ujp. For the simple cells, the plasma 
frequencies are approximately 3 eV, hence these systems 
are highly reflecting for w < 3 eV. 

Eq. [5] is a good starting point for a discussion of the 
trends across Table HIl It basically states that the squared 
plasma frequency is proportional to the product of the 
electron density and the square of the slope of the en- 
ergy bands, both calculated at the Fermi level. In the 
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FIG. 7: Electronic energy band structure of Mgo.75Tio.25H2 
in the simple and super cell along the F-X direction. 



super cells Wp clearly follows the density of states (see, 
Fig. 2]) as a lower concentration of titanium means a lower 
amount of free electrons. In the simple cells this trend is 
not obvious. Subtle changes in the shape of the DOS, i.e. 
bandstructure effects, make for less changes in the elec- 
tron density at the Fermi level. Consistently the plasma 
frequency shows little variation. This, however, does not 
imply that small changes in the slope of the energy bands 
do not play an equally important role here. 

The difference between the simple and super cells of 
the same composition is explained by a decrease in en- 
ergy band slope. The DOS' at the Fermi level are of com- 
parable size, at least for compositions with x — 0.5 and 
X = 0.75. This clearly points to the change in denk/de^ 
as the cause of the significant decrease of cUp. This is 
corroborated by inspection of the band structure, that 
is plotted in Fig. [7] for Mgo.75Tio.25H2 for both the sim- 
ple and super cell. We see that the bandstructure of the 
supercell cannot be understood as a simple back-folding 
of the bands. The randomness in the structure has in- 
duced many interactions between the d-bands, leading to 
a dramatic reduction of their slopes. Going to even larger 
super cells, the effect may become even stronger, and the 
free carrier reflection concomitantly reduced. Such a cal- 
culation goes beyond the present study. It may require a 
different formalism as the reciprocal space description of 
Eq. [5] is bound to break down for truly random systems. 

Interestingly, the flatness of the energy bands also 



makes that the effective mass of the "free" electrons is 
rather high. This goes some way to explain the rather 
high resistivity measured for these systems. The "flat 
bands" also point to a possible localization of carriers. 

For X — 0.875 the reduction of ujp cannot be under- 
stood as a reduction of only 9e„k/9k. From Fig. [4] a 
substantial reduction of the DOS at -EFcrmi is evident. 
Hence reduction of Wp is much stronger than for a; = 0.5 
and X — 0.75. 

To obtain the dielectric functions of the materials we 
just add the interband part of Sect. IIV Al and intraband 
parts obtained from Wp according to Eq. [TJ It was al- 
ready noted above that the impact on the reflection in 
the visible range will be substantial for the simple cells. 
The values of lu^ in the super cells are low enough to only 
induce minor corrections to the interband dielectric func- 
tion. The implications of the corrections will be discussed 
in the next Section. 



V. DISCUSSION 

The calculations on the random alloy super cells clearly 
demonstrate that breaking of the short range order re- 
sults in two important effects: (a) The interband dielec- 
tric function is smoothened and (b) the plasma frequency 
is lowered. The first translates in smoother reflection and 
transmission spectra for the super cells. The second re- 
sults in a lower reflection edge. The simple cells have 
almost full reflection up to 1-2 eV whereas in the super 
cells full reflection only occurs below ^0.3 eV. 

Figure [8] shows the calculated reflection and transmis- 
sion of a film simulating closely the experimental setup as 
used by Borsa et al^ The simulated film consists of 10 nm 
Paladium / 200 nm Mg2;Ti(i_j,)H2 / (fc = 0) quartz sub- 
strate layers, were we use the e(w) calculated in the su- 
per cells for the Mga;Ti(i_^)H2 layer. The plotted range 
is that of a Perkin Elmer Lambda 900 diffraction spec- 
trometer. For comparison; visible light lies in the range 
of 1.65-3.26 eV. The general features and trends as a 
function of Ti content in both reflection and transmis- 
sion are in agreement with the trends observed in the 
experimentsi^ Indeed the reflection is low and our cal- 
culations describe a "black state". With decreasing Ti 
content the reflection decreases and the transmission in- 
creases. 

Since experiments were performed on magne- 
sium/titanium ratios x = 0.7, 0.8, and 0.9, that would 
have required very large super cells, a detailed compar- 
ison can only be made partially. Interpolating between 
the experimental values for 0.7 and 0.8 we can make a 
reasonable comparison with the calculated results for 
X = 0.75. The experimental reflection lies about 5 %Pt. 
lower than the calculated one but the shape of the curve 
is very similar. The experimental transmission drops 
from 6.5%-0% in the range from 0.5-3 eV, it hence lies 
about a factor 3 higher than the calculated one. The 
shapes of the transmission curves, however shows good 
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FIG. 8: (Color online) Reflection and transmission of a 10 nm 
Paladium / 200 nm Mga;Ti(i_^)H2 / {k — 0) quartz substrate 
film. For the dielectric function of Mg2,Ti(i_3;)H2 the results 
for the supercell, including inter and intraband contributions, 
are used. The inset shows experimental results from Ref. |^, 
reproduced with permission, copyright American Physical So- 
ciety (2007). 



agreement. 

Comparing the calculated values for x = 0.875 to the 
experimental results for x = 0.9 we again see somewhat 
lower values for the experimental reflection. Further- 
more, the slight oscillation seen in the calculated spec- 
trum, caused by interference, is much stronger in the ex- 
perimental spectrum. The main difference in the trans- 
mission lies in the energy range below 1 eV where the 
calculated transmission is much larger than the experi- 
mental values. This difference could point to an underes- 
timation of the LUp in the super cell x = 0.875 calculation. 

From the correspondence between the experimental 
and calculated optical spectra we conclude that the 
"black state" can already be explained from moderate- 
size randomized super cells. Put into other words: ran- 
domized models containing only 32 to 64 formula units 
allow for a breaking of the order on a length-scale such 
as to lower the reflection edge and smoothen the spectra. 
This does not necessarily imply that a randomization at 
larger length scales, with a concomitant increase of the 



short range order, e.g. in the coherent crystal picture pro- 
posed by Borsa et al^, would be inconsistent with experi- 
ment. Indeed, the coherent crystal model is supported by 
various observations, e.g. the large positive enthalpy of 
mixing of Ti and Mg. Moreover, for x — 0.875 our model- 
ing seems to underestimate ujp, suggesting that the short 
range disorder may have disrupted the band dispersion 
to much. A full first-principles study of necessarily larger 
models with more short range order is beyond present 
computational capabilities. Such a study is desirable. 
However, with the present randomized models we can 
capture most of the essential physics of Mga;Tii_2;H2. 
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